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Abstract 

The hard X-ray polarimeter POLAR aims to measure the linear polarization of the 50 - 500 keV photons arriving from 
the prompt emission of y-ray bursts (GRBs). The position in the sky of the detected GRBs is needed to determine 
their level of polarization. We present here a method by which, despite of the polarimeter incapability of taking 
images, GRBs can be roughly localized using POLAR alone. For this purpose scalers are attached to the output of 
the 25 multi-anode photomultipliers (MAPMs) that collect the light from the POLAR scintillator target. Each scaler 
measures how many GRB photons produce at least one energy deposition above 50 keV in the corresponding MAPM. 
Simulations show that the relative outputs of the 25 scalers depend on the GRB position. A database of very strong 
GRBs simulated at 10201 positions has been produced. When a GRB is detected, its location is calculated searching 
the minimum of the;^^ obtained in the comparison between the measured scaler pattern and the database. This GRB 
localization technique brings enough accuracy so that the error transmitted to the 100% modulation factor is kept 
below 10% for GRBs with fluence Ftot > 10"^ erg cm"^. The POLAR localization capability will be useful for those 
cases where no other instruments are simultaneously observing the same field of view. 
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1. Introduction 

Gamma-ray bursts (GRBs) are short flashes of y-rays 
likely produced during the creation of a black-hole at 
cosmological distances. In a few seconds a huge amount 
of energy between 10^^ and 10^^ erg is released such 
that GRBs are the most violent explosions in the Uni- 
verse. This prompt emission presents a great variability 
of lightcurves and its spectrum follows a broken power 
law function, often described using the Band model LU- 
The Band model fits the spectrum using four parame- 
ters: an amplitude (A), the low- and high-energy spec- 
tral indexes {a and respectively), and the peak energy 
(Epeak) of the power density spectrum vFy, which repre- 
sents the total energy flux per energy band. 

After the prompt y-ray emission has finished, a long 
lasting afterglow can be observed at various longer 
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wavelengths. The current theoretical picture is that 
GRBs are produced when a massive star collapses at 
the end of its life or when two compact objects merge, 
giving birth in either case to a black hole. The collap- 
sar scenario would be responsible for the long duration 
[2J and the merger scenario for short duration GRBs | 3 1. 
The separation line between short and long GRBs was 
traditionally taken at 2 seconds |4|, but recently sug- 
gested rather at 5 seconds Q. For recent reviews on 
GRBs, see e.g. I^JTl,^. 

Several theoretical models have been proposed to 
explain which processes are responsible for the GRB 
emission. In the standard fireball model | 9 1 the explo- 
sion produces a jet in which expanding shells cause in- 
ternal shocks when a fast shell overtakes a slower mov- 
ing one. In these shocks relativistic electrons are ac- 
celerated, generating through synchrotron emission the 
X- and y-ray photons that we see as the GRB prompt 
emission. The resulting photon polarization is typically 
low (^^ 10%) reflecting the hydrodynamic nature of the 
model. In the electromagnetic model ifTOl the energy 

November 17, 2010 



is extracted from the central engine by an electromag- 
netic field, transported in a plasma to large distances by 
strongly magnetized wind, and dissipated through fast 
magnetic reconnections that generate the GRB radia- 
tion. Photons emitted in this process can present polar- 
ization as high as 50%. In the cannonhall model 1 1 1 1 the 
y-ray emission is generated by inverse Compton scatter- 
ing of soft photons that hit the previously ejected mate- 
rial. The polarization level depends on the opening an- 
gle and Lorentz factor of the jet, predicting therefore all 
possible values between and 100% polarization. 

Precise polarimetry measurements of the GRB 
prompt emission can distinguish between diff'erent theo- 
retical models, off'ering unique information on the emis- 
sion mechanism of the GRBs and on the composition 
and geometric structure of their jets |[l2l[T3l. To date, 
only a few measurements of the prompt GRB polar- 
ization have been performed, all of them with instru- 
ments that had not been designed for this purpose and 
lacked in many cases a good characterization of their 
systematic eff'ects, enough eff'ective area, or good back- 
ground rejection mechanisms for polarimetry. Using 
data from the RHESSI satellite an 80% ± 20% lin- 
ear polarization level from GRB 021206 was reported 
|[T4l . but disproved afterwards (15] [161. By simulat- 
ing the scattering of GRB photons off' the Earth's at- 
mosphere, the polarization levels of two BATSE bursts 
(GRB 930131 and GRB 960924) could be obtained (13 
as >35% and >50% respectively, but the result could 
not be constrained beyond systematics. Also incon- 
clusive were the two measurements (T8| [191 obtained 
using INTEGRAL data from the GRB 041219A. Al- 
though they found a high level of polarization (-60% 
and 98%±33%, respectively) they could not statistically 
claim a polarization detection. 

In view of the power and the lack of precise polariza- 
tion studies of GRBs, several X- and y-ray polarimeters 
have been proposed and are under development. Some 
examples are POLAR |20], GRAPE (Gamma-Ray Po- 
larimeter Experiment, (211) POET (Polarimeters for En- 
ergetic Transients, f22l), CIPHER (Coded Imager and 
Polarimeter for High Energy Radiation, [ 23 1), PHENEX 
(Polarimetry for High ENErgy X rays, 1231), XPOL (25l 
and POLARIX |26|. In addition, polarimeters designed 
for studying fixed sources, like PoGOLite |27|, could 
also measure GRB if they happen to appear in their field 
of view. 

POLAR 1 20 1 is a small and compact instrument de- 
signed to determine the level of linear polarization of 
the 50 - 500 keV photons arriving from the prompt 
emission of GRBs. Totally dedicated to polarimetry, 
the large field of view (-1/3 of the sky) of POLAR 



and its lack of imaging capability would in principle 
force the detector to rely on other instruments to pro- 
vide the location of the observed GRBs. Such a limita- 
tion would reduce the number of GRBs to be measured 
by POLAR depending on whether another GRB detec- 
tor would be observing the same portion of the sky or 
not. To minimize this drawback we have developed a 
method to roughly localize GRBs using only POLAR. 
For this purpose, scalers are attached to the output of the 
multi-anode photomultipliers (MAPMs) that collect the 
light from the POLAR scintillator target. Simulations 
demonstrate that the relative output of those scalers, ac- 
cumulated over the duration of a GRB, shows a depen- 
dency on the source position above POLAR. 

This paper is organized as follows. Section [2] de- 
scribes the POLAR detector and its Monte Carlo sim- 
ulation package. In section [3] we present the working 
principle of the GRB localization method in detail. To 
characterize its capabilities several GRB observations 
have been simulated. The results obtained are presented 
in section [4] where the inffuence on the polarization 
determination of the uncertainty associated to the es- 
timated GRB position is also quantified. Several sys- 
tematic eff'ects that could inffuence the outcome of the 
method are discussed in section[5] including GRB polar- 
ization, background sources, GRB spectral variations, 
and satellite backscattering. Finally, section [6] summa- 
rizes the paper. 

2. Description of POLAR 

POLAR, see ffgure [T] consists of a target of 40 x 
40 plastic scintillator bars (BC400, Saint- Gobain), each 
one 6x6x200 mm^ wrapped in a highly reffective foil 
(Vikuiti, 3M). The optical photons produced when a 
hard X-ray photon interacts with the target are collected 
using 25 MAPMs (H8500, Hamamatsu) whose chan- 
nels are optically coupled to the bottom of the scintilla- 
tor bars. The electrical signal coming from the MAPM 
is then processed by ASIC and FPGA at the front-end 
electronics, and collected by the POLAR central com- 
puter which is responsible for data storage and trans- 
port to ground. The scintillator target is divided into 25 
modules, each consisting of 64 scintillator bars, optical 
coupling, one MAPM, and its corresponding front-end 
electronics, all together enclosed in a thin carbon fiber 
socket. This modular design gives a good mechanical 
stability and facilitates the interchange of modules dur- 
ing the testing phase of the detector. The 25 modules 
are kept together with two aligning frames located at 
the top and at the bottom of the carbon fiber sockets, 
respectively and connected with mechanical fixations. 



The whole target, together with the central computer, 
the power supplies and the rest of the electronics, is 
further enclosed in a box that serves not only as con- 
tainer but also as shield against low energy charged par- 
ticles. Finally, the whole instrument will be mounted 
onto a satellite, to be able to study photons in the en- 
ergy range between 50 and 500 keV, which cannot reach 
the ground because they are absorbed by the Earth at- 
mosphere. A flight opportunity for POLAR on the fu- 
ture Chinese Tian-Gong Space Station is currently un- 
der consideration. An alternative option of flying with 
the International Space Station is evaluated. 




Figure 1 : Scheme of POLAR detector. Left: Exploded view of one 
module from POLAR target. Right: Complete POLAR target, i.e., the 
assembly of 25 modules, with its approximate dimensions. 

An incoming particle is called an event when it de- 
posits more than 5 keV in at least one of the POLAR 
scintillator bars, i.e when it produces a hit. Hard X-ray 
photons arriving from a GRB have a high probability 
of experiencing Compton scattering in the target, gen- 
erating a signal in more than one channel. The POLAR 
trigger logic selects those events with at least two hits. 
When an incoming particle produces more than two hits 
above the 5 keV threshold, the two highest energy de- 
positions are selected oflline. Each of those pairs of hits 
is a selected event in POLAR. The position of the ac- 
tive bars in a selected event is related to the azimuthal 
Compton- scattering angle of the incoming photon. If 



the y-ray emission from the GRB was not polarized, 
the modulation curve, i.e. the azimuthal distribution of 
the ensemble of selected events, is flat. Otherwise it 
follows a sinusoidal curve whose amplitude is the so- 
called modulation factor, and whose phase-shift indi- 
cates the angle of linear polarization of the GRB pho- 
tons ||28l|29l. The modulation factor (yu), divided by the 
response of the detector to a 100% polarized flux (yuioo), 
is the polarization level (EE) of the incoming photons: 
n = jd/jdioo- This principle allows POLAR to mea- 
sure the level of polarization of the prompt emission of 
GRBs. 

The main part of the POLAR front-end electronics 
is devoted to collect the signals produced in the 1600 
MAPM anodes, to implement the trigger, and to trans- 
port the triggered signal to the central computer. The 
trigger condition is fulfilled when the target registers at 
least two hits above the 5 keV single channel thresh- 
old in a time window of ~70 ns. Only the events with 
this characteristic will arrive to the central computer and 
be available for analysis. The probability of registering 
an accidental coincidence is negligible considering the 
counting rates expected during a GRB. 

For the GRB localization method, the electronics 
need to include in addition 25 scalers, one per MAPM, 
counting all photons that pass a higher energy thresh- 
old. If the final electronics design includes a 2-level 
ASIC, it is possible to define a low energy single chan- 
nel threshold of 5 keV to use in the trigger logic, and 
a higher energy threshold of 50 keV per channel, use- 
ful for the GRB localization method implementation. A 
simple "or" operation on the signals that have passed 
the high level threshold, during a time span of about 70 
ns, would be enough to know if at least one hit with 
more than 50 keV energy deposition has been produced 
in the corresponding MAPM. An alternative option for 
the case where a single level ASIC is used would be to 
attach the scaler to the dynode output of each MAPM, 
that provides an analog sum of the 64 channels. All the 
results presented in this paper have been obtained with 
Monte Carlo simulations that reproduced the first ap- 
proach. Nevertheless, independently of which approach 
is taken, the final scaler output for one event is a collec- 
tion of 25 numbers with values or 1 . When accumu- 
lating the scaler output for all events produced during 
a GRB one obtains the number of photons that inter- 
acted above the given high energy threshold inside every 
MAPM. The counting pattern is sensitive to the direc- 
tion of incoming photons since the size of the target cor- 
responds to several absorption lengths. The scalers can- 
not distinguish between GRB photons and background. 
For this reason the influence of the background has been 



studied and is presented in section [5^ The algorithm is 
relatively simple and suitable to be implemented in a 
FPGA or in the small on board processor. The needed 
computing power is negligible as the algorithm has to 
run only infrequently when a candidate GRB has been 
detected. 

An important part of the POLAR design and the eval- 
uation of its polarimetric capabilities is performed using 
GEANT4 |30| simulations. The 1600 scintillator bars 
of the target have been simulated with their standard di- 
mensions and wrapped with a 50 jim thick aluminum 
foil. The target is further placed inside an aluminum 
box of 1 mm thickness that serves as a shield from low- 
energy charged particles and represents the outer enclo- 
sure of the detector. The incoming direction, the spec- 
trum, and the type of particles to be used are generated 
following the user instructions. A specific routine has 
been developed to produce hard X-ray photons follow- 
ing the Band model spectrum d and its parameters. 
The most important physical processes that photons can 
undergo in the detector are taken into account, includ- 
ing polarized Compton scattering, photoelectric eff'ect, 
and pair production, among others. The scintillation and 
optical light collection processes have not been consid- 
ered. The simulation generates a ROOT ISTIl file con- 
taining all the necessary data for the subsequent analy- 
sis: the incoming photon energy, the number and posi- 
tion of the bars fired, the energy deposited at each bar, 
etc. 



3. GRB Localization Method 

The pattern of the scalers output (SO) presented in 
section |2] depends on the position of the GRB above 
the target. A GRB located at POLAR zenith will pro- 
duce an SO with more counts in the entries correspond- 
ing to the central MAPMs. On the other hand, a GRB 
at a large polar angle will present more counts in the 
MAPMs more directly illuminated by the GRB because 
of the exponential interaction probability. Using this re- 
lation it is possible to localize the sources observed by 
POLAR. 

For this purpose a database of SOs has been cre- 
ated using the POLAR Monte Carlo simulation pack- 
age. Non-polarized GRBs have been simulated at nu- 
merous incoming angles. To get enough statistics for the 
comparison with future GRB measurements, we simu- 
lated very strong bursts for the database, with a total 
energy fluence Ftot = 10"^ erg cm"^ and Band spectral 
parameters a = -l.O, JS = -2.5, and Epeak = 200 keV 
m . The selected Band parameters correspond to the ap- 



proximate average values calculated with a selection of 
strong GRBs observed by BATSE |[32ll33]| . 

Defining in spherical coordinates the positions of the 
simulated GRBs in the sky would imply an indetermina- 
tion of the incoming photon azimuthal angle ((py) when 
its polar angle (Oy) approaches 0. To avoid this numeri- 
cal problem the database has been produced in Cartesian 
coordinates: 



X = (l - cos 6y^ • cos (py 
y = (l - cos 6y^ • sincpy 



(1) 



Simulations were uniformly distributed in the (x,y) 
plane, for x and y running between -1 and 1 in steps of 
0.02. We obtain in this way a grid with 10201 nodes. 
Note that, since the detector is symmetric and the in- 
fluence of the satellite behind it is negligible (see sec- 
tion |5.4|), it is enough to simulate one quadrant of the 



grid. The SOs of the other three quadrants can be in- 
ferred from the simulated ones. The selection of the co- 
ordinate system from equation [T] guarantees a uniform 
grid where all points are equally distant from each other 
and the center of the coordinate system corresponds to 
the zenith of the detector. 

The output of the scalers was calculated considering 
only hits with more than 50 keV deposited energy, fol- 
lowing the first approach presented in section [2] Al- 
though there are no intrinsic limitations that prevent the 
detection of weaker energy depositions, the inclusion of 
this threshold is necessary to reduce the spectral depen- 
dence of our localization method (see section |5.3| for 
details). Furthermore, a normalization of the simulated 
SOs is needed due to the dependency of the absolute 
number of counts registered in each MAPM on the to- 
tal flux that the source produced. We scaled the rates 
dividing by the sum of all counts in the given SO. The 
result is a normalized scaler output (NSO) like the ones 
graphically represented in figure [2] The final database is 
constituted by 10201 text files containing the NSO and 
its statistical error for each of the 25 MAPM of POLAR 
target. 

The process to localize the GRBs detected by PO- 
LAR is based on the minimization of the Pearson's 
13411351: 



Mi 



(2) 



where n is the number of MAPMs, c/ is the number of 
counts observed in MAPM /, and m/ is the number of 
counts predicted by a given model for the same MAPM 
/. When POLAR observes a GRB, we calculate c/ as the 
jth gj^^j^y Qf measured SO, and m/ is obtained from the 




Figure 2: Graphical representation of two of the normalized scaler 
output (NSO) patterns that constitute the database. 



database of simulated NSOs. Supposing that the GRB 
was at coordinates (x,y), the predicted number of counts 
in the MAPM is the entry of the database file at 
this position of the sky (m/(x, j)) multiplied by the total 
number of entries (ctot) in the measured SO: 



25 



X' 



(ci-Ctofmi(x,y)) 



(3) 



To localize a GRB we calculate all of 10201 val- 
ues, one at each point of the database grid, obtaining a 
2-dimensional;^^ distribution (see upper part of figurejs] 
for an example done with a burst of same spectral shape 
as the database, but fluence Ftot = 10"^ erg cm"^). The 
minimum of the distribution is found selecting the node 
of the grid with the lowest and fitting two parabo- 
las around it, one in the x-direction and the other in 
the j-direction, using 5 points in each case. The grid 
is fine enough so that the result obtained in this way is 
a good approximation to the absolute minimum of the 
two-dimensional distribution. To perform the parabolic 
fits with ROOT one needs to assign an error to each of 
the points. Since the x^ itself has no associated error 
we assigned a fixed value, equal for all the points. We 
have chosen this value equal to 23, which is the number 
of degrees of freedom of our experiment. On the hor- 
izontal axis the error of each point is the half-width of 
the bin, in our case 0.01. The positioning of the min- 
imum of the parabolas and;^^^(y) corresponds 
to the estimated position of the source. The error on the 
estimation is calculated with the values of the abscissa 
located at;^^^(x)-h 1 and;^^.^(j)-h 1 respectively. Figure 
[3] shows an example of the two-dimensional distribution 
of values, and the two fitted parabolas. 

In figure [s] one can notice that the x^ values do not 
grow monotonously from the minimum, but decrease at 
the four corners of the plane. These four corners of the 
Cartesian space correspond to situations where the GRB 
is below the POLAR detector, i.e. Oy > 90°. The;^^ de- 
creases there because POLAR cannot measure the point 
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Figure 3: Example of minimization. The distribution of in the 
(x,_y) plane is shown in the top, and the parabolic fits of both variables 
in the bottom. The measurement used in this example was produced 
using the POLAR simulation package with a non-polarized GRB of 
Ftot = 10"^ erg cm-2, Qf=-1.0,yS=-2.5, £:peak=200 keV located at PO- 
LAR zenith. 



of interaction of the photons along the length of its bars, 
and is therefore incapable of distinguishing a GRB at 
{6y, (py) from one at (180°-^^, cpy). To avoid this ambi- 
guity, and since GRBs will never appear below POLAR, 
we have introduced an extra condition in the localization 
method that takes it into account. Namely, every time 
that the method gives an output below POLAR it will 
be transformed into the equivalent position above PO- 
LAR and only this one will be given as final result. In 
practice one would actually not need to simulate those 
points of the database. Our reason to do it is the numer- 
ical advantage of working with a uniform and square 
grid that can be treated as a regular matrix. 

4. Method verification 

Simulations were made to produce example measure- 
ments with which we could test the localization method. 
The GRBs simulated as examples were non-polarized 
standard strong GRBs. Hereafter, standard strong 
stands for a GRB with total energy fluence Ftot = 10"^ 
erg cm"^ and Band spectral parameters |1| a = -LO, 
/3 = -2.5, and Epeak = 200 keV, from which POLAR 
will detect about 22000 photons in the whole duration 
of the burst. Note that the spectal shape, i.e. the Band 
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Figure 4: Results from simulated non-polarized standard strong GRBs. Top left: Linearity plot in x, i.e. reconstructed x vs. input x for several 
fixed _y values. Each point on this graph represents the average of the 10 simulations made at each position. Middle left: Spread of the reconstructed 
X, i.e. the standard deviation of the 10 reconstructed x values. Bottom left: Bias in the reconstruction, i.e. the difference between the reconstructed 
average x value and the input x. The plots on the right are the equivalent to the left ones, but calculated for y at several fixed x positions. In all the 
plots lines of different styles have been drawn to guide the eye. 



parameters, is the same as used for the database bursts, 
but the total fluence is here ten times smaller. We ex- 
pect, considering the BATSE catalog and the POLAR 
field of view, to observe about 12 such standard strong 
GRBs per year |36|. To facilitate the visualization we 
present only GRB positions above POLAR and in the 
quadrant of the space where x > and j > 0. The be- 
havior on the other three quadrants is equivalent due to 
the symmetry of the detector. 



We compared the x and y values taken as input in 
the measurements with the ones reconstructed using 
the minimization method. Ten samples were run 
at each selected {x,y) position to study the repeatabil- 
ity and the spread of the reconstructed results around 
each input value. Figure [4] shows the performance of 
the GRB localization method. The localization proce- 
dure works well reconstructing both x and y without ap- 
parent bias within a root-mean- square (r.m.s.) smaller 



6 



than 0.07. The average minimum obtained for the 
380 measurements simulated to produce figure]?] was 
= 22.9 ± 6.6. Taking into account that there 
are 23 degrees of freedom, this corresponds to a reduced 
= 0.996 ± 0.287. When performing the equiv- 
alent GRB localization in spherical coordinates Qy could 
be determined within a r.m.s ^ 4° and 0^ within a r.m.s 

5° for Qy > 30°. As comparison, the localization ac- 
curacy of the Fermi-GRB burst monitor (GBM) ITTl is 
for bright bursts < 1.5° using a detailed offline analysis, 
and the Swift Burst Alert Telescope (BAT) \3W\ provides 
positions to an accuracy of 1' - 4' within 20 seconds of 
GRB observation. 

When a GRB is detected by POLAR its level of po- 
larization (n) is determined from the amplitude (ji) of 
the modulation curve and the 100% modulation factor 
(yUioo) using: 11 = yu/yUioo (see section |2]). The error on 
the polarization level can be calculated as: 

where ctu, cT/j, and ct^^qq are the errors on H, jj. and yUioo, 
respectively. 

The second term in equation [4] depends only on the 
strength and level of polarization of the observed GRB. 
Since yUioo varies with the burst position and spectral 
shape, n must be calculated using a yUioo derived from 
simulations of a GRB with the same spectrum and po- 
sition in the sky as the observed one. The highest yUioo 
is obtained when the GRB is located at POLAR zenith 
and it diminishes when photons come from POLAR side 
(see figure [5]). When the location of the GRB is un- 
known, the value of yuioo needed to calculate U is also 
unknown. Therefore the uncertainty on the position of 
the GRB is translated to an uncertainty in the value of 
jJiQQ and, through the first term of equation |4j further 
into an error on the resultant level of polarization. 

One can calculate the error of yUioo with standard error 
propagation: 

For each of the positions simulated as examples, the 
spread and bias of the reconstructed values were deter- 
mined. When the spread of the reconstructed coordi- 
nates was larger than their bias, the spread was intro- 
duced as (Tx and o-y in equation [5] Otherwise, the value 
of the bias was taken. To calculate the partial deriva- 
tives of the yUioo with respect to x and y we fitted the 




Figure 5: Distribution of 100% modulation factor (//loo) vs. x and y, 
calculated for a standard strong GRB (see text). The empty area at the 
right of the image corresponds to positions where the GRB would be 
below POLAR. Such a situation will never happen during flight and 
therefore was not considered in the analysis. 

2-dimensional distribution from figure |5] with the fol- 
lowing polynomial function: 

f^iooi^, y) = ao -\- aix -\- a2X^ -\- a^y + a/^y^ + asxy{6) 

whose partial derivatives are: 

1^ = ai+ 2a2X -h asy 
(7) 
^ = a3y + 2a^y -h a^x 

From equations [7] and |5] one can calculate the error 
induced on yUioo by the uncertainty in the location of the 
GRB. The final result can be seen in figurej6]and repre- 
sents the value of the first term in equationpFcr^^QQ //^loo) 
at various source locations. The error caused to yUioo is 
always smaller than 6% of its absolute value. 




Figure 6: Error caused to //loo due to the uncertainty on the localiza- 
tion of GRBs, divided by //loo- 

The total error in the measured polarization is given 
by equation [4] Its second term (cr^/yu) depends on the 
strength and level of polarization of the GRB. In figurejT] 
we present an example calculated for standard strong 



bursts of fluence Ftot = 10"^ erg cm"^, located at PO- 
LAR zenith, and with various levels of polarization. For 
this calculation we have substituted cr^ioo/A^ioo = 0-10 
in equation |4j instead of the 0.06 previously mentioned, 
because 10% was the maximal error found on yUioo for 
such a standard strong GRB when variations of the 
spectral shape were taken into account (see section [53] ). 
One can notice that the error due to the position uncer- 
tainty is negligible for low levels of polarization, with 
respect to the error associated to the modulation fac- 
tor itself. At large levels of polarization the localiza- 
tion contribution to the error is more important. The 
final error in 11 is in that case cru < 15% of the abso- 
lute polarization level. The values represented in fig- 
ure [7] should be considered for bursts of total fluence 
Ftot = 10"^ erg cm"^ as upper limits to the error on the 
polarization level. 
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Figure 7: Upper limits to the error on the polarization level for a GRB 
of total fluence Ftot = 10~^ erg cm~^. These plots have been calcu- 
lated using the upper limit of the error inflicted to //loo due to the lo- 
calization procedure. Left: Relative error of the modulation factor for 
difl'erent levels of polarization {dotted), together with the final relative 
error on IT (solid), where o-^^^Iii\qq = 0.1 has been included. The in- 
crease due to the localization uncertainty is visible at high polarization 
values. Right: Calculated polarization level {points), together with the 
simulation input values {line), plotted as a reference. For polarization 
levels above 60% the uncertainty in the GRB position is the largest 
source of error. 



5. Systematic effects 

The results presented in figure |4] represent the best 
case example, since the measurements to which the 
method was applied were produced in the same con- 
ditions as the localization database, and considering 
POLAR as an instrument with a perfectly uniform re- 
sponse. During a real observation there are several is- 
sues that could influence the outcome of the GRB local- 
ization method. We will discuss here the most important 
ones: asymmetries produced by GRB polarization, dif- 
fuse photon background, GRB spectral variations, satel- 
lite backscattering, fluctuations in light-collection efl&- 
ciency, and MAPM non-uniform sensitivity. 



5.7. GRB polarization 

The database of the localization method was con- 
structed simulating unpolarized sources. Although the 
polarization introduces an asymmetry in the POLAR 
hit pattern, it does not strongly afl'ect the output of the 
scaler used for the localization procedure. Some ex- 
amples were made with simulations of standard strong 
GRBs, 100% polarized in the direction parallel to PO- 
LAR X-axis, to confirm that statement. As with the 
non-polarized examples, 10 simulations were produced 
at each position of the sky and the spread and bias in the 
reconstruction of their locations were calculated. It was 
clearly seen that the outcome of the localization proce- 
dure was not aff'ected by the polarization of the source. 
Both the spread of the reconstructed values and the bias 
of their average were within the errors equal to the non- 
polarized case. 

5.2. Diffuse y-ray background 

Outside of the South Atlantic Anomaly the largest 
source of background that will aff'ect POLAR on orbit is 
the diff'use y-ray photon background 1201 . This source 
of background illuminates approximately isotropically 
the POLAR instrument. Since the surface of the most 
external MAPMs is more exposed, the SO produced 
by the y-ray diff'use photon background presents higher 
number of counts for the external MAPMs and lower in 
the internal ones. When analyzing the SO of the total 
signal, i.e. GRB-hbackground, the result of the localiza- 
tion method may be diff'erent than the obtained from the 
GRB alone. POLAR will store data from before, dur- 
ing, and after the GRB, making possible the subtraction 
of the diff'use y-ray background. We will discuss here 
how accurately the interpolation and subtraction need to 
be made so that the localization method is not aff'ected. 

The spectrum of the y-ray diff'use photon background 
has been parametrized 1391 using the function: 

3-60 keV : 

f(E) = 7.877£-0-29^-^/4i-i3 

> 60 keV : 

/(£) = 0.0259(1) + 
+0.504(1)-'" + 
+0.0288(1) 

with f(E) in units of , ^ • An increment 

keV cm s sr 

of 10% in the absolute normalization factor has been 
needed to fit recent measurements of the cosmic hard 
X-ray spectrum performed by INTEGRAL instruments 

mni. 



To study the background influence it is important to 
consider the duration of the burst. UnHke GRBs, where 
for a constant fluence the flux diminishes when the du- 
ration increases, the difl'use background produces a flux 
approximately constant with time. Since we accumu- 
late the output of the scaler during the whole GRB, 
the shape of its lightcurve has no influence on the lo- 
calization method. For simplicity we considered both 
the lightcurves of the burst and the background as flat. 
Figure [8] shows the lightcurve and spectra of the difl'use 
photon background together with a short GRB of 1 sec- 
ond duration and standard Band parameters located at 
POLAR zenith. The equivalent plots for a 20 seconds 
long GRB are presented in figure [9] 
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was subtracted from the previous data successively with 
an excess of 8%, 10%, 20%, and 30% over its average 
value. The whole procedure was performed twice, using 
a 1 second and 20 seconds duration GRB, to reproduce 
a short and long GRB, respectively. 

We observed that, while for a short GRB the back- 
ground did not influence the performance of the local- 
ization method, its contribution was very important in 
the case of the long GRB, inducing an error of -50% 
in yuioo if not subtracted. If the background is overes- 
timated by a 8% of its real level, the error transmitted 
to yUioo was cr^^QQ <10%, increasing to ct^^qq <14% for a 
10% overestimation, and to ct^^qq ~20% and (T^^qq -30% 
for a 20% and 30% overestimation, respectively. 

As a conclusion, background subtraction is important 
not only for the polarimetric performance of POLAR, 
but also for its capability of localizing GRBs. In the 
case of short GRBs the difl'use y-ray background can 
be neglected, but for 20 seconds long GRBs the back- 
ground calculation must not difl'er from the real back- 
ground level by more than 8%. When the background 
is not strongly variable, it is sufficient to average two 
background measurements one taken immediately be- 
fore and the other immediately after the GRB. 



Figure 8: Comparison between the flux of a short (1 sec) standard 
strong GRB and the difl'use photon background. Left: Assumed 
Hghtcurve of the GRB signal (solid) and the difl'use photon back- 
ground (dotted). All photons producing a hit above 5 keV in POLAR 
target have been included. Right: Spectrum for the same two sources. 
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Figure 9: Comparison between the flux of a long (20 sec) standard 
strong GRB and the difl'use photon background. Left: Assumed 
lightcurve of the GRB signal (solid) and the difl'use photon back- 
ground (dotted). All photons producing a hit above 5 keV in POLAR 
target have been included. Right: Spectrum for the same two sources. 

We have simulated the diffuse y-ray background us- 
ing equation [8] with a 10% higher normalization factor 
as suggested in ll40ll , and we calculated its SO. The re- 
sult was added to the SOs of a series of standard strong 
GRBs as the ones presented in figure |4] Then, a sec- 
ond diff'use y-ray background was simulated and its SO 



5.3. Variations of the GRB spectral shape 

The GRBs used for the NSO database have been sim- 
ulated taking the Band parameters corresponding to a 
standard strong GRB, summarized in the first row of ta- 
ble[T] The NSOs and SOs depend on the spectrum since 
the penetration length of photons in the POLAR target 
is an exponential function of their energy. To study the 
possible influence of difl'erent GRB spectra on the local- 
ization method, several series of simulations have been 
performed. We selected from the BATSE catalog |33 1 
the bursts with energy fluence Ftot ~ 10"^ erg cm"^ 
which had the Band parameters most difl'erent from our 
standard ones. The selection of GRBs studied is pre- 
sented in table [U 

We made a series of simulations locating each 
BATSE burst at several positions in the sky and applied 
to them the localization method to study the spread and 
bias of the reconstructed values, in the same way as pre- 
sented in section |4] When using a very low 5 keV en- 
ergy threshold to the scalers, it was observed that the 
modification of the GRB spectral parameters introduced 
in some cases a large bias on the reconstructed x and y 
coordinates. The sign of the bias changed from GRB 
to GRB, and its absolute value was the largest (taking 
values up to 0.25) when a and Epeak strongly diff'ered 
from the values used for the database, while variations 
on did not aff'ect it significantly. The main reason is 



GRB 


Ftot (xlO-^) 


a 




Epeak 




[ erg cm"^] 






[keV] 


Standard 


1.0 


-1.0 


-2.5 


200 


911104 


1.1 


-0.6 


-1.95 


265 


920718 


1.1 


-1.09 


-3.06 


192 


930922 


2.0 


-1.40 


-2.88 


94 


951016 


1.4 


-1.61 


-2.02 


118 


980828 


1.7 


-0.25 


-2.08 


223 


981130 


1.1 


-0.54 


-2.25 


649 



Table 1: Details of the spectrum of the standard strong GRB com- 
pared to the bursts from BATSE catalog selected for studying the in- 
fluence of the spectral shape into the localization procedure. From the 
BATSE catalog we selected the six bursts that, having a similar to- 
tal energy fluence, presented the most diff'erent Band parameters with 
respect to our standard GRB. 



that most of the photons that POLAR detects, and which 
mostly define the shape of the SO, have low energies. 
Low energy photons do not penetrate far inside POLAR 
target, producing hits only on the side where the GRB 
is coming from. Although this makes them in principle 
very good to localize a GRB, the fact that POLAR can- 
not measure the spectral shape precisely enough makes 
it impossible to correct for the bias in the reconstructed 
coordinates for GRB spectrally very diff'erent from our 
standard one. 

One possible solution would be to measure the 
GRB spectrum with a small spectrometer that could be 
mounted somewhere close to POLAR. In such a case 
it would be enough to produce a new NSO database 
with the observed burst spectrum and apply the local- 
ization method to determine the GRB position with a 
few degrees uncertainty. When no spectrometer is avail- 
able, the only alternative to mitigate the problem is, 
as we have done, to raise the energy threshold of the 
scaler, which in principle could be as low as 5 keV, to 
50 keV. In this way we largely reduce the sensitivity of 
the method to spectral variations, so that the bias of the 
reconstructed coordinates is below 0.15 and the error 
inflicted into yUioo is below 10% for all GRBs from ta- 
ble [T] We show in figure [TO] the case of GRB 981130, 
for which the method performance was the worst of the 
six examples (the absolute values of the coordinates bias 
were the largest). 

Unfortunately, when reducing the method sensitiv- 
ity to spectral variations, we cut out a large number of 
counts from the SOs and NSOs, diminishing the statis- 
tical power. As a consequence, the error produced on 
yUioo from the uncertainty in the localization of GRBs 
with fluence Ftot 



parameters as the database bursts, can reach 14%. For 
weaker GRBs our technique will not be able to provide 
the position of the source with sufficient accuracy, espe- 
cially taking into account the spectral variability. 

5.4. Satellite backs cattering 

Photons coming from the GRBs can reach POLAR 
indirectly after being backscattered in the spacecraft 
where the polarimeter is mounted. We have simulated 
a very simple description of the Chinese Tiang-Gong 
Space Lab as a "bottle-like" assembly of two cylinders 
(1.65 m radius, 4 m high, and 1.4 m radius, 3.2 m high, 
respectively) joint by a 1.4 m long trunked cone. All 
volumes were simulated in aluminum with a thickness 
of 4 cm, so that the total mass of the spacecraft is ^^8 
tons. POLAR was placed outside the Space Lab at 
12.5 cm from its trunked cone part, pointing with the 
front of its target to the zenith. The GRB photons were 
uniformly illuminating an area of 2 m radius with PO- 
LAR at its center. A series of simulations with standard 
strong GRBs at several positions in the sky, as done for 
figure]?! were produced and the localization method was 
applied. The results were found to be within errors the 
same as when the satellite was not considered. The 50 
keV energy threshold applied in the scalers is rejecting 
most of the photons coming to POLAR from the satel- 
lite, since they have lost energy in the backscattering 
process. 

5.5. Statistical fluctuations in light collection 

Plastic scintillators are a good choice to perform hard 
X-ray polarimetry because the probability that a photon 
of this energy range experiences Compton scattering in 
the target is very high. The price one pays in exchange 
is a low energy sensitivity, not enough to make spectro- 
scopic studies. GEANT4 simulations of optical photon 
tracking show that about ^^60% of the optical photons 
produced in a POLAR scintillator bar never reach the 
MAPM The detector will be calibrated taking this 
into account but, due to the statistical fluctuations in the 
number of optical photons collected, the spectral reso- 
lution will be poor. We have simulated the scintillation 
process in a single bar and fired it with bunches of pho- 
tons at various energies. For each run a gaussian fit of 
the photo-peak spectral line was performed. We found 
an expression to relate the cr of this gaussian fit with the 
incoming energy of the photon, i.e., the broadening (cr) 
of the lines as a function of the deposited energy (Edep, 
expressed in keV): 



= 10 ^ erg cm ^ and same spectral 



10 



cr = 0.215 • Edep[keV] -h 2.953 keV. 



(9) 
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Figure 10: Results from GRB 981130, simulated non-polarized and without background. Top left: Linearity plot in x, i.e. reconstructed x vs. input 
X for several fixed _y values. Each point on this graph represents the average of the 10 simulations made at each position. Middle left: Spread of 
the reconstructed x, i.e. the standard deviation of the 10 reconstructed x values. Bottom left: Bias in the reconstruction, i.e. the difference between 
the reconstructed average x value and the input x. The plots on the right are the equivalent to the left ones, but calculated for y at several fixed x 
positions. In all the plots lines of different styles have been drawn to guide the eye. 



The main consequence of this Hne-broadening for PO- 
LAR measurements is that the signal collected at the 
MAPM anode does not correspond exactly to the de- 
posited one. Therefore, the 50 keV threshold introduced 
at the scalers will be blurred. 

The influence of the spectral line broadening into the 
GRB localization method has been studied. Standard 
strong GRBs were simulated at difl'erent positions in the 
sky as it was done to produce figure]?] Before applying 



the localization method to these GRBs, each value of 
the energy deposition Edep was substituted by a number 
that had been randomly generated following a Gaussian 
distribution with mean Edep and sigma the result from 
equation [9] The application of the localization method 
to the modified data gave the same results as when Edep 
was perfectly determined, with no increment on the co- 
ordinates uncertainty. 



5.6. Sensitivity variations between MAPM anodes 

Two photons of the same energy can produce differ- 
ent signals if they come to two MAPM channels that 
do not have the same sensitivity. According to the 
H8500 MAPM data-sheet |42l sensitivity differences 
up to a factor of 3 can exist between the anodes of 
the same MAPM. We have observed that the sensitivity 
usually varies following a smooth function with a max- 
imum close to one of the comers of the MAPM, and 
monotonously decreases when moving away from that 
point. We have assigned to each of the 64 elements of a 
MAPM a number in the range from 0.4 to 1, distributed 
according to an example data-sheet, but with variations 
of the order of 5% at each channel. Such a sensitivity 
mask was produced for each of the 25 MAPMs of PO- 
LAR, taking random orientations so that the corner with 
maximum sensitivity was not always at the same side of 
the target. In this study we have assumed that the av- 
erage sensitivity is the same for all 25 MAPMs. The 
40 X 40 mask created in this way was applied as a mul- 
tiplicative factor to the energy depositions of a series of 
standard strong GRBs like the ones from figure |4] No 
change in the localization uncertainty with respect to the 
results shown in figure [4] was found. 

6. Summary and conclusions 

Relevant information on the magnetic composi- 
tion, geometric structures, and emission mechanism of 
GRBs, not accessible by spectroscopic and photometric 
approaches, can be obtained through polarization mea- 
surements of the ORB prompt emission. POLAR aims 
to determine the level of linear polarization of GRBs 
in the energy range from 50 to 500 keV by measuring 
the azimuthal distribution of the photons that Compton- 
scatter in its plastic scintillator target. To keep the de- 
sign as simple and compact as possible, POLAR is to- 
tally devoted to polarimetry and is not optimized for 
imaging and spectroscopic purposes. The level of polar- 
ization of a GRB can only be determined when knowing 
its position in the sky. 

We have presented a method, using only POLAR, 
able to estimate the position of GRBs with enough pre- 
cision so that the error transmitted to the 100% modu- 
lation factor due to the location uncertainty is kept be- 
low 10% for strong GRBs of total fluence Ftot > 10"^ 
erg cm"^. A scaler is connected to the output of each 
MAPM to record the number of GRB photons that pro- 
duce at least one hit with energy above 50 keV. The out- 
puts of the 25 scalers, accumulated over the whole du- 
ration of the GRB, present a pattern that depends on the 
direction where the photons are coming from. 



We have performed several series of simulations 
to estimate the capabilities of the GRB localization 
method and the influence of its uncertainty on the mea- 
sured level of polarization. When simulating strong 
non-polarized GRBs with a fluence Ftot = 10"^ erg cm"^ 
and a similar spectral shape to the ones used to create 
the database, the position coordinates could be deter- 
mined with a maximal error of 0.07, equivalent to being 
able to determine the position in spherical coordinates 
(Oy, (py) within ~ 5° for Oy > 30°. This would inflict 
an error in the yUioo of the measured GRB of less than 
6%. Considering GRBs of similar strength but differ- 
ent spectral shapes, the coordinates could be determined 
with a maximal error of 0.14, keeping the error on yUioo 
below 10%. For weaker GRBs of fluence Ftot ~ 10"^ 
erg cm"^ the error in yUioo will be below about 20%. 

Several effects that could potentially modify the per- 
formance of the GRB localization technique have been 
studied. It was found that the influence of the diffuse 
photon background will be negligible for short GRBs, 
but it needs to be subtracted from long ones. Its esti- 
mation in the later case should be done with a preci- 
sion better than 8% around the real value. Neither the 
level of polarization of the source, nor the presence of 
a spacecraft behind POLAR that backscatters part of 
the GRB photons show any influence in the output of 
the localization method. The limited energy resolution 
due to light collection inefficiencies and MAPM non- 
uniformities will blur the 50 keV energy threshold but 
this will not diminish the localization capabilities of our 
method. 

The GRB localization technique will provide enough 
accuracy to allow for the measurement of the GRBs po- 
larization level with only a small increase on its error. 
For a GRB of total fluence Ftot = 10"^ erg cm"^, the 
added error is negligible at low polarization levels, com- 
pared with the one from the measured modulation factor 
itself. At large polarization levels the error from the lo- 
calization is larger than the one from the measured mod- 
ulation factor, but all together reaching no more than 
15% of the measured polarization. 

The method presented here can be experimentally 
tested when uniformly illuminating the complete PO- 
LAR target with a radioactive source to be placed at 
various locations. Such a kind of test will be performed 
when the POLAR engineering qualification model will 
be ready. In flight, if other instrument can provide the 
position of one or several GRBs detected by POLAR, 
a direct comparison with the output of the localization 
method can be done. In general, the precision in the lo- 
calization provided by instruments with imaging capa- 
bilities will be much better than that of POLAR alone. 



and preferred to it. The POLAR localization capabil- 
ity will be useful for those cases where no other instru- 
ments are simultaneously observing the same field of 
view. 
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